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Abstract. We study the distribution of maxima (Extreme Value Statistics) for sequences 
of observables computed along orbits generated by random transformations. The underly- 
ing, deterministic, dynamical system can be regular or chaotic. In the former case, we will 
show that by perturbing rational or irrational rotations with additive noise, an extreme 
value law will appear, regardless of the intensity of the noise, while unperturbed rotations 
do not admit such limiting distributions. In the case of deterministic chaotic dynamics, 
we will consider observables specially designed to study the recurrence properties in the 
neighbourhood of periodic points. The exponential limiting law for the distribution of 
maxima is therefore modified by the presence of the extremal index, a positive parameter 
not larger than one, whose inverse gives the average size of the clusters of extreme events. 
The theory predicts that such a parameter is unitary when the system is perturbed ran- 
domly. We perform sophisticated numerical tests to assess how strong is the impact of 
noise level, when finite time series are considered. We find agreement with the asymptotic 
theoretical results but also non-trivial behaviour in the finite range. Finally, we focus on 
a prototypical dissipative system such as the Henon map with standard values for the 
parameters. In a previous paper, we showed that it is possible to relate the parameters 
describing the type of the limiting extreme value laws with the local dimension of the 
attractor of the system. Here, we find that while it is well known that adding an arbitrary 
weak noise results into smoothing the invariant measure, it is far from trivial to detect 
such an effect. Our results suggest that in many applications where finite datasets can be 
produced or analysed one must be careful in assuming that the smoothing nature of noise 
prevails over the underlying deterministic dynamics. 
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1. Introduction 

The main purpose of this paper is to study the extremal behaviour of randomly perturbed 
dynamical systems. By extremal behaviour we mean its statistical performance regarding 
the existence of Extreme Value Laws (EVLs), or in other words, the existence of distribu- 
tional limits for the partial maxima of stochastic processes arising from such systems. In 
many aspects of natural and social sciences and engineering, the statistical properties of 
the extremes of a system are usually relevant tied with actual risk assessment. This is an 
element of why the theory of extremes has received such a great deal of attention over the 
years. The motivation for considering randomly perturbed systems follows from the fact 
that, very often, dynamical systems are used to model human activities or natural phenom- 
ena and the fact that errors made by observations usually take a random character which 
can be well described by the random perturbation formalism. On the other hand, random 
noise is often added in numerical modelling in order to represent the lack of knowledge on 
(or practical impossibility of representing) some of the processes taking place in the sys- 
tem of interest, often characterised by small spatial and/or temporal scales, whose explicit 
representation is virtually impossible. See [31] for a comprehensive discussion on this issue 
in a geophysical setting. Also, noise is often considered as a "good" component to add 
when performing numerical simulations, because of its ability to smoothen the invariant 
measure and basically remove unphysical solutions [28] . 

Very recently, in [2], the authors proved the first results (up to our knowledge) regarding 
the existence of EVLs for piecewise expanding systems which are randomly perturbed by 
additive noise. They observed that adding noise has a "smoothing" effect in terms of elim- 
inating all possible clustering of exceedances. The distributional properties of the maxima 
of stationary sequences is driven by the appearance of exceedances of high thresholds. For 
independent and identically distributed (i.i.d.) processes, the exceedances appear scattered 
through the time line. For dependent processes this may not necessarily hold and clustering 
of exceedances may occur. The Extremal Index (EI), which we will denote by < 9 < 1, is 
a parameter that quantifies the amount of clustering. No clustering means that 9 = 1 and 
strong clustering means that 9 is close to 0. For deterministic hyperbolic systems, it has 
been shown, in [TTJ [121 (HI 12], that a dichotomy holds: either an EI < 9 < 1 is present 
at the (repelling) periodic points, or the EI is always 1 at every other point. We remark 
that Hirata [IB] had already observed the two types of behaviour but not the dichotomy. 
Coming back to the "smoothing" effect of adding noise, to be more precise, in [2], the 
authors show that, for the randomly perturbed systems considered, the EI is always 1. 
The main tool used there was decay of correlations against L 1 observables, which held for 
the unperturbed and perturbed systems considered. We remark that this property implies 
an excellent mixing behaviour of both systems. 

One of the main achievements of this paper is the extension of the results in [2], by 
showing that adding noise still has a smoothing effect when the unperturbed system is not 
mixing at all, or even worse: when the system is actually periodic. To be more precise, 
for unperturbed systems, we consider rotations on the circle (irrational or not) and show 
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that by adding absolutely continuous noise with respect to the Haar measure, we can 
always prove the existence of EVLs with the EI being 1 everywhere. The analysis of these 
systems will be made using Fourier transforms to compute decay of correlations of specific 
observables. Observe that for the original unperturbed systems the lack of mixing does 
not allow to prove EVLs in the usual sense. Hence, here, we have a more drastic transition 
which motivates the question of whether it is possible to distinguish numerically the real 
nature of the underlying (unperturbed) systems when we look at the extremal statistics of 
the randomly perturbed data. 

Another main goal of this paper is to support the analytical discussion with numerical 
simulations devised to show how, for finite samples, the extremal behaviour may or not 
follow the asymptotic results. This has obvious relevance in terms of applications. One 
of the key aspects surfacing from the numerical analysis of extremes that we carried out 
is that the estimation of some dynamical and geometrical properties of the underlying 
physical measure strongly depends on a combination of the intensity of the noise and the 
length of the data sample. 

There are two main experimental issues deriving from stochastic perturbations. The first 
concerns with the general claim that the addition of noise has a "smoothing" effect over the 
physical measure and enhances the chaoticity of the underlying deterministic dynamical 
systems. As we have mentioned, we expect to see this in terms of extremal behaviour, at 
least asymptotically. This claim will be weighed up against the simulations showing that 
we may need to integrate the system for a very long time to observe any changes in the 
measure structure. Moreover, for higher dimensional systems featuring stable and unstable 
manifolds, the stochastic perturbations should be directed on the stable direction so that 
the "smoothing" effect of the measure is more effective. 

The second issue is related to the level of noise needed to modify the deterministic prop- 
erties and, in this case in particular, the sensitivity of the extremal type behaviour to the 
the noise level. 

To discuss this latter issue we will use the fact that the numerical round off is comparable 
to a random noise on the last precision digit [20]. This observation will allow us to claim 
that, for systems featuring periodic or quasi periodic motions (as the rotations on the 
circle), such a noise level is not sufficient for producing detectable changes regarding the 
observed extremal behaviour. 

On the other hand, the analysis for generic points of chaotic systems, such as the full shift 
on three symbols, will show that the intrinsic chaotic behaviour of the map dims the effect 
of the noise and makes the perturbed system indistinguishable from a deterministic one. 
If, instead, the analysis is carried out at periodic points, the effect of the perturbation 
will be easily noticeable in the disruption of periodicity and its consequent clustering of 
exceedances, which can be detected by estimation of the EI. 

The disruption of periodicity will be further analysed through simulations on the quadratic 
map. We will show, for finite sample behaviour, how the EI depends on the intensity of 
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the noise. We will conclude the analysis suggesting a way to analyse the disruption of 
periodicity in general systems where the nature of the noise may be unknown. 

We will also present some numerical simulations on the Henon map, which features a sta- 
ble and an unstable directions, homoclinic tangencies and the coexistence of two different 
attractors. In a previous paper we have shown that it is possible to relate the EVL pa- 
rameters to the local dimension of the attractor of a system possessing a singular invariant 
physical measure [29]. This latter experiment will be the gateway to go beyond the maps 
presented in this paper and to suggest a general procedure to analyse issues related to 
the stochastic perturbations of dynamical systems in order to frame the analysis carried 
out here in a more general setup. In fact, we find that it is far from trivial to detect the 
smoothing effect of noise on the invariant measure when finite time series are considered. 

We would like to remark that such considerations have great relevance in the context of 
Axiom A systems and having statistical mechanical applications in mind when addressing 
the problem of the applicability of the fluctuation-dissipation theorem (FDT). While in the 
deterministic case the non-smooth nature of the measure along the stable manifold makes 
it impossible to apply straightforwardly the FDT (see discussion in [331 Ell ESI 123, ED] ) , the 
addition of even a very small amount of random forcing in principle "cures" the singularities 
of the measure and makes sure that in principle the FDT can be used to related forced 
and free motions (see [U [21]). Our results suggest that one must be careful in assuming 
that such a desirable effect is really detectable when finite datasets are analysed. 



Our main purpose is to study the extremal behaviour of randomly perturbed dynamical 
systems. This will be done by analysing the partial maxima of stationary stochastic pro- 
cesses arising from such systems. We will start by presenting the main concepts regarding 
the extreme values on a general framework of stationary stochastic processes. Then we 
will turn to the construction of stationary stochastic processes deriving from dynamical 
systems and their random perturbations by additive noise, which will be the object of our 
study. 

2.1. Extreme values — definitions and concepts. Let Xq,Xi, ... be a stationary sto- 
chastic process. We denote by F the cumulative distribution function (d.f.) of X , i.e., 
F(x) = P(Xo < x). Given any d.f. F, let F = 1 — F and uf denote the right endpoint of 
the d.f. F, i.e., up = sup{x : F(x) < 1}. We say we have an exceedance of the threshold 
u < up whenever 



2. Extreme values, the system and the perturbation 



U(u) := {X > u} 



(2.1) 



occurs. 



We define a new sequence of random variables M\ , M2 , . . . given by 

M n = max{A , . . . , X n _i}. 



(2.2) 
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Definition 1. We say that we have an Extreme Value Law (EVL) for M n if there is a 
non-degenerate d.f. H : K — > [0, 1] with -f^(O) = and, for every r > 0, there exists a 
sequence of levels u n = u n {r), n = 1, 2, . . ., such that 

nP(X > u Tl ) — y t, as n — » oo, (2.3) 

and for which the following holds: 

P(M n < u n ) -> H(t), as n ->■ oo. (2.4) 

where the convergence is meant in the continuity points of H(r). 



The motivation for using a normalising sequence u n satisfying (2.3) comes from the case 
when X ,Xi, . . . are independent and identically distributed (i.i.d.). In this i.i.d. setting, 
it is clear that P(M n < u) = {F{u)) n , where F is the d.f. of X , i.e., F{x) : = P(X < x). 



Hence, condition (2.3) implies that 



P(M n < un) = (1 - P(X > u n )) n ~ (l - ^ e" r , (2.5) 

as n — 7- oo. Moreover, the reciprocal is also true (see [2U Theorem 1.5.1] for more details). 
Note that in this case H(r) = 1 — e~ T is the standard exponential d.f. 

When Xo, Xi, X2, ■ ■ ■ are not independent but satisfy some mixing condition D(u n ) intro- 
duced by Leadbetter in [22J then something can still be said about H. Let i 7 ^,...,^ denote 
the joint d.f. of X h ,. . . ,X in , and set F iu ,„ >in (u) = F ht _ >in (u, ...,«). 

Condition (D(u n )). We say that D(u n ) holds for the sequence Xq,Xx, ... if for any inte- 
gers %x < ■ ■ ■ < i p and ji < . . . < jh for which j\ — i p > m, and any large n e N, 

\F iu ...,i p , h ,..., jk {u n ) - F h ^ ip (u n )F jlr ..j k (u n )\ < 7(n,m), 

where 7(71, m n ) )■ 0, for some sequence m n = o(n). 



If D(u n ) holds for X ,X!, . . . and the limit (2.4) exists for some r > then there exists 



< 6 < 1 such that H{r) = e~ 6r for all r > (see [231 Theorem 2.2] or [21 Theorem 3.7.1]). 

Definition 2. We say that Xq,X\, . . . has an Extremal Index (EI) < 9 < 1 if we have 
an EVL for M n with H(r) = e~ 9r for all r > 0. 

The notion of the EI was latent in the work of Loynes [26] but was established formally by 
Leadbetter in [23]. It gives a measure of the strength of the dependence of X , Xx, . . ., so 
that 6 = 1 indicates that the process has practically no memory while 8 = 0, conversely, 
reveals extremely long memory. Another way of looking at the EI is that it gives some 
indication on how much exceedances of high levels have a tendency to "cluster" . Namely, 
for 9 > this interpretation of the EI is that 9~ l is the mean number of exceedances of a 
high level in a cluster of large observations, i.e., is the "mean size of the clusters". 

In the dependent shown in [22], if D(u n ) holds and in addition an anti clustering 

condition D'(u n ) (see definition below) also holds, one can show that the EI is 1, which 
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means that H(r) = e _T , as in the independent case. However, since the rates of mixing 
for dynamical systems are usually given by decay of correlations of observables in certain 
given classes of functions, it turns out that condition D{u n ) is too strong to be checked for 
chaotic systems whose mixing rates are known only through decay of correlations. For that 
reason, motivated by Collet's work [8], in [H] the authors suggested a condition _D 2 (wn) 
which together with D'(u n ) was enough to prove the existence of an exponential EVL 
(H(t) = e~ T ) for maxima around non-periodic points z. In fact, [HI Theorem 1] states 
that if the following conditions hold for X ,X ly . . . then there exists an EVL for M n and 
H(t) = 1 - e- T . 

Condition (D 2 (u n )). We say that D 2 (u n ) holds for the sequence X ,X 1 , ... if for all £, t 
and n 

\P(X >u n f] max{X t , . . . ,X t+e ^ < u n }) - P(X > u n )F(M e < u n )\ < 7 (n,t), (2.6) 

where 7(^,t) is decreasing in t for each n and n^(n,t n ) — > when n — > oo for some 
sequence t n = o{n). 

Now, let (fcn)neN be a sequence of integers such that 

k n — > oo and k n t n = o{n). (2.7) 
Condition (D'(u n )). We say that D'(u n ) holds for the sequence Xq, X%, X 2 , ... if there 



exists a sequence {A;„}„ e N satisfying (2.7) and such that 

[n/k n \ 

lim n V P(X > u n ,X j > u n ) = 0. (2.8) 

n— >oo ' ^ 

Condition D 2 (u n ) is much weaker than the original D(u n ), and it is easy to show that it 
follows easily from sufficiently fast decay of correlations (see [HI Section 2]). 

In the results mentioned above, condition D'(u n ) prevented the existence of clusters of 
exceedances, which implies for example that the EVL was a standard exponential H (r) = 
e~ r . However, when D'(u n ) does not hold, clustering of exceedances is responsible for the 
appearance of a parameter < 9 < 1 in the EVL which now is written as S(t) = e~ 9r . 

In [17] , the authors established a connection between the existence of an EI less than 
1 and periodic behaviour. This was later generalised for rare events point processes in 
[TB"] . Under the presence of periodic phenomena, the inherent rapid recurrence creates 
clusters of exceedances which makes it easy to check that condition D'{u n ) fails (see [T7J 
Section 2.1]). This was a serious obstacle for the application to dynamical systems since the 
theory developed up to [TTJ was based on Collet's important observation that D'{u n ) could 
be used not only in the usual way as in Leadbetter's approach, but also to compensate 
the weakening of the original D(u n ), which allowed the application to chaotic systems 
with sufficiently fast decay of correlations. To overcome this difficulty, in [IT] , the authors 
considered the annulus 

Q p (u) := {X >u, X p < u} (2.9) 
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resulting from removing from U(u) the points that were doomed to return after p steps, 
which form the smaller ball U(u) H f~ p (U(u)). In [TTJ [TB], the occurrence of Q p (u) was 
named an escape since it corresponds to the realisations escaping the influence of the 
underlying periodic phenomena and exit the ball U(u) after p iterates. Then, the main 
crucial observation in [T7j is that the limit law corresponding to no entrances up to time 
n into the ball U(u n ) was equal to the limit law corresponding to no entrances into the 
annulus Q p {u n ) up to time n (see [TTJ Proposition 1]). This meant that, roughly speaking, 
the role played by the balls U(u) could be replaced by that of the annuli Q p (u), with the 
advantage that points in Q p (u) were no longer destined to return after just p steps. 

Based in this last observation and motivated by the behaviour at repelling periodic points, 
in [17], the authors proposed certain general conditions (for general stationary stochastic 
processes) designed to prove the existence of an EI less than 1 under the presence of some 
sort of underlying periodic phenomenon. The first one establishes exactly the type of 
periodic behaviour assumed, namely: 

Condition (SP Pi g(u n )). We say that X , Xi, X 2 , . . . satisfies condition SP p ^(u n ) for p G N 
and 9 e [0, 1] if ' 

lim sup F(Xj > u n \X > u n ) = and lim F(X P > u n \X > u n ) (1 - 0) (2.10) 

n— >oo X<j<p n— >oo 

and moreover 

lim P(X > u n , X p > u n , X 2p > u n , . . . , X ip > u n ) = 0. ( 2 -H) 



n— >oo 

i=0 



Condition (2.10), when 9 < 1, imposes some sort of periodicity of period p among the 
exceedances of high levels u n , since if at some point the process exceeds the high level 
u n , then, regardless of how high u n is, there is always a strictly positive probability of 
another exceedance occurring at the (finite) time p. In fact, if the process is generated by 



a deterministic dynamical system / : X — t X and / is continuous then (2.10) implies that 
z is a periodic point of period p, i.e., f p (z) = z. 

The next two conditions concern to the dependence structure of X ,Xi, . . . and can be 
described as being obtained from D 2 (u n ) and D'(u n ) by replacing balls by annuli. Let 
&n(u n ) := f]j=o f~^(Q( u n)°)- Note that while the occurrence of the event {M n < u n } 
means that no entrance in the ball {X > u n } has occurred up to time n, the occurrence 
of £? n (u n ) means that no entrance in the annulus Q(u n ) has occurred up to time n. 

Condition (D p (u n )). We say that D p (u n ) holds for the sequence Xq, X\,X 2 , ... if for any 
integers i, t and n 



|P(Q P ,oM n ^W) - P(Q p , (« n ))P(J p ,o/K))| < 7 (M 



where j(n,t) is non increasing in t for each n and n / y(n,t n ) — > as n — > oo for some 
sequence t n = o(n). 
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This condition requires some sort of mixing by demanding that an escape at time is an 
event which gets more and more independent from an event corresponding to no escapes 
during some period, as the time gap between these two events gets larger and larger. The 
main advantage of this condition when compared to Leadbetter's D(u n ) (or others of the 
same sort) is that it follows directly from sufficiently fast decay of correlations as observed 
in Section 3.3 of [IT] , on the contrary to D{u n ). 

Assuming D p {u n ) holds let {k n ) n& n be a sequence of integers such that 

k n — > oo and k n t n = o{n). (2-12) 
Condition (D' (u n )). We say that D' Ju n ) holds for the sequence X ,X 1 ,X 2 , ... if there 



exists a sequence {k n } n( z N satisfying (2.12) and such that 



[n/k n ] 

lim n V P(Qp,oK) H Q p j(u n )) = 0. (2.13) 



This last condition is very similar to Leadbetter's D'(u n ) from [23], except that instead of 
preventing the clustering of exceedances it prevents the clustering of escapes by requiring 
that they should appear scattered fairly evenly through the time interval from to n — 1. 

In [TTl Theorem 1], it was proved that if a stationary stochastic process satisfies conditions 
SP p> e(u n ), D p {u n ) and D' (u n ) then we have an EVL for M n with H(r) = e~ dT . 



2.2. Stochastic processes arising from randomly perturbed systems. To simplify 
the exposition we will consider one dimensional maps / on the circle S 1 or on the unit 
interval I (from now on we will use I to identify both spaces), provided, in the latter case, 
that the image f(I) is strictly included into /. We will perturb them with additive noise, 
namely, we introduce the family of maps 

f u (x) = f(x) +u 

where we have to take the mod-1 operation if the map is considered on the circle. The 
quantity u> is chosen on the interval Q £ = [—e,e] with distribution i? e , which we will take 
equivalent to Lebesgue on Q £ . Consider now an i.i.d. sequence Uf., k e N taking values 
on the interval Q e and distributed according to d £ . We construct the random orbit (or 
random transformation) by the concatenation 

fZ(x) = f^of LUn _ 1 o...of UJl (x). (2.14) 

The role of the invariant measure is now played by the stationary measure fi £ which is 
defined as 

J J <l>{fJ,x)) dfJL e (x) d$®((j) = J <f)(x) dfj, E (x), 
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for every <\> G L°° to be intended with respect to the Lebesgue measure Leb)[j The 
previous equality could also be written as 



U £ <p dfi £ 



where the operator U £ : L°° — > L°°, is defined as (U £ (fi)(x) = L <fi(fuj(x))d$ £ and it is called 
the random evolution operator. 

The decay of correlations of the perturbed systems could be formulated in terms of the 
random evolution operator. More precisely, we will take two non-zero observables ,_0 and 
ip and we will suppose that <fi is of bounded variation with norm 
the correlation integral is 

1 



\bv an d ip G LA Then 



Cox e ((j),ip,n) : = 



U £ (f) ipdfi £ - / <p dfi £ / tp d\x 



\BV\\W\\Ll 

In the following we will use the product measure P := \i e x # 



We are now in condition of defining the time series Xq, X\, X2,. . . arising from our system 
simply by evaluating a given observable (p : I — > R U {+00} along the random orbits of the 
system: 

X n = V o£, for each n£l, (2.15) 

We assume that ip achieves a global maximum at z G J; for every u < ip(z) but sufficiently 
close to <p(z), the event {y G / : (p(y) > u} = {X > u] corresponds to a topological 
ball "centred" at z and, for every sequence (u n ) n( z^ such that u n —> <p(z), as n — > 00, the 
sequence of balls {U n } n( z^ given by 

U n := {X > u n } (2.16) 
is a nested sequence of sets such that 

H U n = {4- (2-17) 

nSN 

As explained in [TJ] the type of asymptotic distribution obtained depends on the chosen 
observables. In the simulations, we will cover the observables such that 

¥>(•) = <7(dist(.,s)), (2.18) 
where dist denotes a certain metric chosen on / and g is of one of the following three types: 

(1) gi(y) = — log(y) to study the convergence to the Gumbel law. 

(2) g 2 {y) = y~ 1/a a G R a > for the Frechet law. 

(3) ^ 3 (V) = C - y x l a a G R a > C G R for the Weibull distribution. 



This choice is dictated by the fact that the stationary measure will be equivalent to Lebesgue in all 
the examples considered below 

2 0f course we could do other choices, but this two spaces will play a major role in the subsequent theory. 
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The sequences of real numbers u n = u n (r), n = 1,2, . . ., are usually taken to be as one 
parameter linear families like u n — y/a n + b n , where j/Gl and a n > 0, for all nGl In 
fact, in the classical theory, one considers the convergence of probabilities of the form 

F(a n (M n -b n ) <y). (2.19) 

In this case, the Extremal Types Theorem says that, whenever the variables X t are i.i.d, if 
for some constants a n > 0, b n , we have 

F(a n (M n -b n ) <x) -+G{x), (2.20) 

where the convergence occurs at continuity points of G, and G is non degenerate, then G 
belongs to one of three extreme values types (see below). Observe that r depends on y 
through u n and, in fact, in the i.i.d. case, depending on the tail of the marginal d.f. F, we 
have that r = r(y) is of one of the following three types (for some a > 0): 

n(y) = e-» for y e R, r 2 (y) = y~ a for y > and r 3 (y) = {-y) a for y < 0. (2.21) 

In [211 Theorem 1.6.2], it were given sufficient and necessary conditions on the tail of the 
d.f. F in order to obtain the respective domain of attraction for maxima. Besides, in 
[2U Corollary 1.6.3] one can find specific formulas for the normalising constants a n and b n 
so that the respective extreme limit laws apply. We used these formulas to perform the 
numerical computations in this paper. 

Remark 1. We emphasise that, for i.i.d. sequences of random variables, the limiting dis- 
tribution type of the partial maxima is completely determined by the tail of the d.f. F. 
For the stationary stochastic processes considered here, if an EI 9 > applies, then the 
same can still be said about the limiting distribution type of the partial maxima: namely, 
it is completely determined by the tail of the d.f. F. This statement follows from the 



equivalence between (2.3) and (2.5) and the definition of the EI. However, we also quote 
[2~it Corollary 3.7.3] because we will refer to it later: 

If for X , Xi, ... we have an EI 6 > 0, then if we considered an i.i.d. sequence Z , Zx, . . . 
so that the d.f. of Z is F, the same as that of X , and let M n = max{Z , . . . , Z n _i}, then 
the existence of normalizing sequences (a n ) ng N and (6 n ) nS N for which 

lim P(a n (M„ - bn) <y)= G{y) 

n—>-oo 

implies that 

lim F(a n (M n -b n )<y) = G e (y), 

n— >oo 

and the reciprocal is also true. Moreover, since by [2U Corollary 1.3.2] G e is of the same 
type of G, we can actually make a linear adjustment to the normalizing sequences (a n )„ g N 
and (6 n )„ e N so that the second limit is also G. 

Remark 2. From Remark [TJ in order to determine the type of extremal di strib ution G 
(recall that G(y) = e T ^ y \ where r(y) is of one of the three types described in (2.21)) which 



applies to our stochastic processes X ,X 1 , . . ., one needs to analyse the tail of the d.f. F. 



The choice of the observables in (2.18) implies that the shape of g determines the type of 



EXTREME VALUE STATISTICS FOR DYNAMICAL SYSTEMS WITH NOISE 



11 



extremal distribution we get. In particular, for observables of type gi we get an extremal 
law of type e T % for i — 1, 2, 3. (See [13 Remark 1] for more details on this correspondence). 
While the type of the extremal distribution is essentially determined by the shape of the 
observable, in the cases when types 2 and 3 apply, i.e., the Frechet and Weibull families of 
distributions, respectively, the exponent a is also influenced by other quantities such as the 
EI and the local dimension of the stationary (invariant) measure /i e (/i). In particular, when 
such measure is absolutely continuous with respect to Lebesgue and its Radon-Nikodym 
derivative has a singularity at z, then the order of the singularity also influences the value 
of a. 



3. Numerical experiments 



For the numerical computations we will use the approach already described in [9]. It 
consists in considering unnormalised maxima selected using the block Maxima approach. 
Once computed the orbit of the dynamical systems, the series of the gi observables are 
divided into m blocks of length n observations. Maxima of the observables g^ are selected for 
each blocks and fitted to a single family of generalised distribution called GEV distribution 
with d.f.: 



Fg(x; v, a, k) = exp 



1 + K 



x — V 



a 



-1/k' 



(3.1) 



which holds for 1 + k(x — u)/a > 0, using v 6 R (location parameter), cr > (scale 
parameter) and k G R is the shape parameter also called the tail index: when k — > 0, 
the distribution corresponds to a Gumbel type ( Type 1 distribution). When the index is 
positive, it corresponds to a Frechet (Type 2 distribution); when the index is negative, it 
corresponds to a Weibull (Type 3 distribution). In this fitting procedure, as we have seen 
in [9], the following relations hold: 

for type 1, n = 0, v = b n , cr = l/a n ; (3.2) 
for type 2, k = 1/a, v = b n + l/a n , a = n/a n ; (3.3) 
for type 3, k = — 1/a, v = b n — l/a n , a = —n/a n . (3.4) 

For the discrete maps like f^x), we select a value for z and repeat the following operations 
for various values $e of the noise intensity: 

(1) We construct an empirical physical measure of the system by performing a long 
run and saving the obtained long trajectory. The realisations / £Jn _ 1 ,-- - , f Ul are 
constructed by taking the with the uniform distribution in Q e . 

(2) We select 500 initial conditions according to the physical measure described above; 
such initial conditions will be used as initial conditions (x variables) to generate 
500 realisations of the stochastic process. 
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(3) For each of the 500 realisations, we obtain an orbit containing r = m ■ n = 10 6 , 10 7 
data. 

(4) We split the series in m = 1000 bins each containing n = 1000 (or n = 10000) ob- 
servations. These values are chosen in agreement with the investigations previously 
carried out in [9], since they provide enough statistics to observe EVL for maps 
which satisfy the D2 and D' conditions. 

(5) We fit the GEV distribution for the g>j, i — 1, 2, 3 observables and the maxima and 
minima and study the behaviour of selected parameters. 

In the forthcoming discussion we will present the results for the shape parameters n(gi) 
and for the location parameter cr(gi). This choice is convenient, as explained in [29j ITU] , 
since these parameters do not depend on the number of observations n and therefore their 
asymptotic expected value does not change when r is modified. For the one dimensional 
maps considered we expect to obtain the following asymptotic parameters for systems 
satisfying the mixing conditions described in the former sections (see [TO]): 

K ( 9l ) = K(g 2 ) = l/a K,(g 3 ) = -l/a a{g 1 ) = 1 (3.5) 

The inference procedure follows [9] where the authors have used a Maximum Likelihood 
Estimation (MLE) that is well defined when the underlying physical measure is absolutely 
continuous. In the applications described in [291 ITU] an inference via a L-moments proce- 
dure has been preferred as it provides reliable values for the parameters even when the d.f. 
corresponds to fractal or multi fractal measures. However, since the L-moments procedure 
does not give information about the reliability of the fit^J hereby we exploit the MLE pro- 
cedure since it helps to highlight the situations where the fitting does not succeed because 
of the poor data sampling. 

The following numerical analysis is constructed in a strict correspondence with the theo- 
retical setup: the role of the simulations is to follow step-by-step the proofs of the theorems 
trying to reproduce the analytical results. Whenever it occurs, the failure of the numerical 
approach is analysed in relation to the length of the sample and the intensity of the noise, 
just to highlight which parameters are crucial for the convergence procedure. 

We focus on simple discrete maps as they already contain interesting features regarding the 
convergence issues. We start by analysing the rotations on the 1-D torus where the condi- 
tions D2{u n ) and D'(u n ) are not satisfied, but they could be restored under perturbation. 
In the case of expanding maps of the interval, for which the existence of EVL has been 
already shown both analytically and numerically, we show that the computation of the EI 
reveals to be crucial to discriminate whether the system was randomly perturbed or not. 



The L-moments inference procedure does not provide any confidence intervals unless these are derived 
with a bootstrap procedure which is also dependent on the data sample size [B]. The MLE, on the other 
side, allows for easily compute the confidence intervals with analytical formulas [32 . 
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Eventually, the analysis of the Henon map, shows how the presence of another basin of 
attraction (infinity in this case) may affect the EVL regardless the initial condition chosen. 

Besides the convergence issues, one of the main problems that we face regarding the esti- 
mation of the EI during the realisation of the numerical simulations is the fact that the 
fitting procedure used in [9], "hides" the EI. This is the case because, as we mentioned 
in Remark [TJ when quoting [2U Corollary 3.7.3], the normalizing constants a n and b n can 
be chosen (after a simple linear rescaling) so that we get exactly the same limit law G 
of the corresponding i.i.d. process instead of G e . This prevents the detection of an EI 
< 6 < 1. Basically, our original procedure simply selects the best fitting constants from 
the sample of maxima available, which means that the EI does not surface at all. In order 



to overcome this difficulty, instead of applying a blind fitting, we use relations (3.2)-(3.4) 



(more specifically (3.2)), the knowledge about the stationary measure of the particular 
systems we considered and the information provided by [T5J (TTJ |2] to compute a priori 
the normalizing sequences so that we can capture the value of the EI. This adjustment to 
the procedure revealed also very effective and the results met completely the theoretical 
predictions. 



4. Rotations on the circle 



In this section we will study the rotations of the circle and show that the effect of adding 
noise is enough to create enough randomness in order to make EVLs appear when they do 
not hold for the original system. This will be the content of Theorem [T] below. Because 
this nice statistical behaviour appears solely on account of the noise, it is not surprising 
that the numerical simulations show that for very small noise one needs a large amount of 
data to the detect the EVLs. 



4.1. Analytical results. In [2J, Theorem C], the following result has been proved essen- 
tially for piecewise expanding maps randomly perturbed like above. 

To be more precise, let us suppose that the unperturbed map / is continuous on the circle 
and that 

(i) the correlation integral Cor e (0, ip,n) decays at least as rT 2 \ 

(ii) u n satisfies: nP(X Q > u n ) — > r, as n — > oo; 

(iii) U n = {X > u n } verifies f] neN U n = {z}; 

(iv) There exists r\ > such that d(f(x),f(y)) < i] d(x,y), where d(-, •) denotes some 
metric on /, 

then the process Xq,X 1} ■ ■ ■ satisfies D 2 (u n ) and D'(u n ), and this implies that the EVL 
holds for M n so that H{r) = e~ T . The proof strongly relies on the decay of correlations for 
L 1 (Leb) observables. 



The result is even more general and applies to multidimensional maps too, but for our concerns, 
especially for rotations, the 1-D case is enough. We will discuss later on about generalisation to piecewise 
continuous maps 



14 D. FARANDA, J.M. FREITAS, V. LUCARINI, G. TURCHETTI, AND SANDRO VAIENTI 



In this section we show that such a result holds also for rotations (irrational or not), 
perturbed with additive noise. The key observation is that in the proof of [21 Theorem C], 
the observables entering the correlation integral are characteristic functions of intervals 
(see also [21 Remark 3.1]) and for such observables it is possible to prove an exponential 
decay of correlations for perturbed rotations by using the Fourier series technique. 

In what follows we identify S 1 = K/Z. 

Theorem 1. Let f : S 1 — > S 1 be a rotation of angle a G IR, i.e., f(x) = x + a mod 1. 
We perturb f additively so that the random evolution evolution of an initial state x G S 1 
is given by (2.14), with $ £ denoting the uniform distribution on [— e,e]. Let X , X±, . . . be 



a stationary stochastic process generated by the random evolution of such f , as in \2.15 
Let (u„)neN be a sequence such that items (ii) and (Hi) above hold. Then, the process 
X Q ,Xi, . . . satisfies conditions D 2 (u n ) and D'{u n ) which implies that there exists an EVL 
for M n with EI equal to 1, i.e., H{t) = e~ T . 

For ease of exposition, we will change slightly the notation for the random perturbation of 
/. We write them in this way: 

f £ ^(x) = x + a + e(, modi (4.1) 

where £ is a random variable uniformly distributed over the interval [—1,1] and therefore 
of zero meanjf] Let us observe that 



f 3 ei (x) = x + ja + e(£ 1 + --- + Z j ) 



(4.2) 



In this case, it is straightforward to check (just by using the definition), that the stationary 
measure coincides in this case with the Lebesgue measure, Leb, on [0, 1] and it is therefore 
independent of e. 

Let us first establish that the correlation of the right functions decays exponentially fast 
under the random evolution of the system. 

Lemma 4.1. Under the assumptions of Theorem [71 if = xa and ip = \b, where x 
denotes the characteristic function, B = U e l=1 Bi, for some £ G N and A, B\, . . . , Bi C S 1 
are connected intervals, then 



C 



U^(i/j)(f)dx — / ipdx I (j)dx 
n Jo Jo 



< 4e" J ' e lo s( 27r ) 



(4.3) 



as long as e 2 < 1 — log2/log(27r). 



Proof. We begin by writing the modulus of the correlation integral, C, e , as follows: 



C 



3,£ 



dx- 



2-i 



i • • • 



i f 1 f 1 

d£,jip(f^.Ax))4> ~ \ 4>dx / <\>dx 
i s Jo Jo 



5 With respect to the previous notations, we changed J\ into [— 1, 1], u = e£, with £ 6 [ — 1, 1] and finally 
$ E becomes d£ over [—1,1]. 
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2-nikx 



Note that given A = [a, b] and B\ = [ai, bj\, I — 1, ...,£, we have 



= / Xa e 
'o 



-2nikx 



dx 



^—2-irika ^—2irikb 



2nik 



and also 



Z" 1 V € e" 

Jo 

which implies that 

< 1/1*1- 

Plugging the Fourier series of <f>,if) into Cj )6 , we obtain 



2-Kikai Q—2mkbi 



2nik 



(4.4) 



C 



3,B 



-17 i f m — 1 — 



fcez/{o} 



22 ^k4>-kS 3 {ke) 

k&/{0} 



sin(27nr) 



where S(xj 
and \S(x)\ < -j^, \x\ > 1. 
We now continue to estimate Cj e as 



. For this quantity we have the bounds: 15(2;) | < e x lo &( 27r ), | x | < ^ 



l/e oo 
-fc 2 e 2 jlog(27r) 



C ie <2j]^Ui0- fc |e- 



k=l 



+ E m 

fc=l/e 



(2vrfe)^ 



(4.5) 



Using (4.4), we have 



C- < 2e- £2jlog(2 " ) V — + 2 V — 

J£ " ^ fc2 + ( 27r )i 2^ fc 2 

fc=l V 7 fc=l/e 

Using the inequality ^2 k=A p < ;f — |>, we finally have 

C7 ij6 < 2(2 - £ ) e -^^(2n) + 4ee -jlog(2^) < 4e -i £ 2 log(2^) 

which is surely verified when e~ je2 lo §( 27r ) > 2e" j log ( 27r \ namely e 2 < 1 — log 2/ log(27r). □ 

Proof of Theorem^ We now check the conditions D2(u n ) and D'(u n ) using the previous 
estimate for the exponential decay of random correlations for characteristic functions. We 
will use a few results established on the afore mentioned paper [2]. We begin with D 2 (u n ) 
and we consider, with no loss of generality, the observable g(x) = — log(d(x, z)), where 
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d(-, ■) is the distance on the circle and z is a given point on the circle. The event U n := 
{g > u n } = {X > u n } will be the ball B e -u n (z). We now introduce the observables 

<fi( x ) = X{X >u n } = X{g>u n }, 
4>(X)=J X {g , go f e(l ,..., g of^<u n } d t 1 - 

It has been shown, in [21 Sect. 4.1], that proving condition D2(u n ) can be reduced to 
estimating the following correlation 

J [l £ (g> U n , g O f eg < Un , . . . , Q O f £ f' X < U n ) - 

/ > «.) df JM9< «., 9 °U<u„...,go /« < „„) J? 



Since all the maps f s g are globally invertible and linear, the events (g > u n ) and (g o < 
u n , ••• ,g° < u n ) are a connected interval and a finite union of connected intervals, 



respectively. Therefore we can apply directly formula (4.3) with an exponential decay 



in t so that any sequence t n = n K , with < k < 1, will allow to verify the condition 
n^{n, t n ) —> 0, when n — > oo. 

The computation of D'{u n ) is a bit more lengthy. Let us fix a sequence ct„ — > oo so that 
a n = o(k n ). Next, we introduce the quantity R-(A) as the first return of the set A into 
itself under the realisation £. As in [21 Sect. 4.1, just before eq. (4.5)], we have: 

[n/k n \ \n/k n \ 

n F (Un n 4 J (t/„)) < n X] ^ef/ n , e U n ,Rl(U n ) > a n }) 

j = l «n 



[n/fc n J 



+ n ^ P({(z,£) : a; G U n ,Rl(U n ) < a n }) 



(4.6) 



i=i 



It can be shown that the measure of all the realisation such that R^{U n ) < a n is bounded 
by a constant times Leb(U n ) times a n (we use here the fact that 7] in item (iv) above is 1 



in our case). Therefore the second term in (4.6) is bounded by a constant times: 



n 
k n 



Leb(U n ) 2 a n - 



This reduces to r 2 ^, whose limit is 0, as n — > oo, given our choices for a n and k n . We 

fan 

now consider the first term on the right hand side of the inequality (4.6). We can obtain 
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the following result: 

\n/k n \ 



J2 P({(z,0 -x e U n , f[(x) e U n ,Rl(U n ) > a n }) < 



n 



n 



U™ J 

Leb(U n ) 2 + n \ C je 



]=a n 



(4.7) 

where the correlation integral is computed with respect to the characteristic function of 
U n . Next we use the decomposition (4.5) for Cj e ; we also observe that the product of the 
modulus of the fc-Fourier coefficient of \u n f° r itself gives 

' (l-cos27rJfeLeb(l7 n )). 



2n 2 k 2 



Using the notation \U n \ = Leb(C/ n ), we obtain (note that we discard the first term in (4.7) 
since it goes to zero because n 2 \U n \ 2 ~ r 2 ): 



L^J l/e 

(4.7) < n { 2 E - ™s2nk\U n \)e 



k 2 e 2 j log(27r) 



j=a„ k=l 



°° 1 1 

+ 2£^(l-™2rf|«,|)p SiF }. 

l/e V ' 



We consider the first term into the bracket: an upper bound can be defined as follows: 



j=a n k=l x 



n 



1 ( \ — cos27rk\Ur, 



When n — > oo, n\U n \ ~ r and the inner series goes to zero. The second piece in (4.7) could 
be bounded as before by 



Lfen J 2 00 Y 

Un \ (2^)7 

j=a„ K ' l/e 



1 — cos27r£;|[y ?) 
\LL\ 



which converges to zero as well. 



□ 



4.2. Numerical Investigation. In this section we describe the results obtained from the 
numerical investigation of the stochastically perturbed rotation map following the proce- 
dure describe in Sect. 3. The results are displayed in FigjTJwhere the green lines correspond 
to experiments where we have chosen a bin length n = 10 4 , whereas the blue lines refer 
to experiments where the chosen bin length is n = 10 3 . The red lines indicate the values 
of the parameters predicted by the theory for a 1-D map satisfying the mixing conditions 
described above. We set e = 1CT P to analyse the role of the perturbations on scales ranging 
from values smaller than those typical for the numerical noise to 0(1) [11]. 
The solid lines display the values obtained by averaging over the 500 realisations of the 
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stochastic process, while the error bars indicate the standard deviation of the sample. Fi- 
nally, with the dotted lines we indicated the experiments where less than 70% of the 500 
realisations produce a statistics of extremes such that the empirical d.f. passes success- 
fully the non-parametric Kolmogorov-Smirnov test [25J when compared to the best GEV fit. 

Even though the theory guarantees the existence of EVL for the rotations perturbed with 
an arbitrarily weak noise (e.g. when compared with the numerical noise), the simulations 
clearly show that EVLs are obtained when considering small but finite noise amplitudes 
only when very long trajectories are considered. The quality of the fit improves when larger 
bins are considered (compare blue and green lines in Figjl]). This is in agreement with the 
idea that we should get EVL for infinitely small noises in the limit of infinitely long sam- 
ples. In our case, EVLs are obtained only for e > 10 -4 , which is still considerably larger 
than the noise introduced by round-off resulting from double precision, as the round-off 
procedure is equivalent to the addition to the exact map of a random noise of order 10~ 7 

[201 H]. 

This suggests that is relatively hard to get rid of the properties of the underlying determin- 
istic dynamics just by adding some noise of unspecified strength and considering generically 
long time series: the emergence of the smoothing due to the stochastic perturbations is 
indeed non-trivial when considering very local properties of the invariant measure as we 
do here. It is interesting to expand the numerical investigation discussed here in order to 
find empirical laws connecting the strength of the noise perturbations with the length of 
the time series needed to observe EVLs. 



5. Expanding Maps 

In [2] it is proved that expanding maps on the circle admit EVLs when they are perturbed 
with additive noise (see, in the just quoted reference, Corollary 4.1 for smooth maps and 
Proposition 4.2 for discontinuous maps on I). The proof is completely different of the one 
we produced for rotations in Theorem [T} The proof there relies on the fact that expand- 
ing maps have exponential decay of correlations for BV-functions against L 1 observables. 
Moreover, as we already anticipated, the proof is not limited to continuous expanding 
maps, but it also holds for piecewise expanding maps with a finite number of discontinu- 
ities, provided that the map is topologically mixing. It is interesting to point out that 



under the assumptions (items (i)-(iv)) in the beginning of Subsection 4.1, we could prove 
convergence to the e~ T law with the observables gi previously introduced independently of 
the choice of the point z. We should remind that for deterministic systems and whenever z 
is a periodic point of prime period p, then the limit law is not e~ T any more, but rather e~ er , 
where < 9 < 1 is the extremal index discussed above. Therefore, obtaining experimen- 
tally an extremal index with unitary value for an observable ip(-) computed with respect 
to a periodic point z definitely points to the fact that the system is stochastically perturbed. 
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If a non-periodic point z is considered, it seems relevant to investigate whether stochastic 
and deterministic systems exhibit any sort of differences when looking in detail into the 
resulting EVLs. We consider the following map: 

fet(x) = 3x + e£ modi (5.1) 

perturbed with additive noise as in the rotation case discussed above. The stationary mea- 
sure for such a map is the uniform Lebesgue measure on the unit interval independently 
of the value of e. The numerical setup is the same as in the case of the rotations: we 
begin by taking a non-periodic point as the center of our target set; the results are shown 
in picture [2] It is clear that the stochastic perturbations do not introduce any changes in 
the type of statistical behaviour observed for a non-periodic point z and no differences are 
encountered even if the number of observations in each bin is increased. This is compatible 
with the idea that the intrinsic chaoticity of the map overcomes the effect of the stochastic 
perturbations. Summarising, extremes do not help us in this case to distinguish the effect 
of intrinsic chaos and the effect of adding noise. 

It seems more promising the problem of the form of the EVLs for periodic points. As 
discussed above, we cannot use the usual fitting procedure for the GEV, since it always 
renormalises in such a way that the extremal index seems to be one (see discussion above). 
Instead, in order to observe extremal indices different from one, we have to fit the series of 
minimum distances to the exponential distribution by normalising a priori the data. The 
inference of the extremal index 9 (here already the extremal index) has been obtained via 
MLE. The normalisation applied consists only in multiplying the distances by the factor 
2n as we deal with a constant density measure. In Fig. [3] we present the results for the 
extremal index obtained taking the periodic point z = 1/2 of prime period 1 for which 
6 = 2/3 |^| The results clearly show that we are able to recognise the perturbed dynamics 
as the extremal index goes to one when e increases. Interestingly, the rupture with the 
values expected in the purely deterministic case is observed only for relatively intense 
noise. Moreover, when longer time series are considered (green experiment), the stochastic 
nature of the map becomes evident also for weaker perturbations. Finally, it is clear that 
the numerical noise (corresponding to e ~ 10~ 7 ) is definitely not sufficiently strong for 
biasing the statistics of the deterministic system. 

6. Quadratic map 

Following the analysis on the extremal index carried out in the previous section, we study 
numerically the quadratic map 

fe^{x) = 1 - ax 2 + e£ 



We recall that 6 = 6{z) = 1 — \detD(f p (z)\, where z is a periodic point of prime period p, see [T71 
Theorem 3]. 
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in the two cases a = 0.014 and a = 0.314. For these values of a, the deterministic dynamics 
is governed by the existence of an attracting fixed point (whose value depends on a); for 
both choices of a the fixed point is taken as z in the experiments. The stationary measure 
of such perturbed quadratic maps is absolutely continuous with respect to Lebesgue [3]. 
Extremes of this map have been already studied in the unperturbed case in 1 13 . but 
for chaotic regime (a belonged to the Benedicks-Carleson set of parameters). Here, we 
proceed exactly as described in the previous section, by fitting the series of minimum 
distances to an exponential distribution after normalizing the data using the factor 2n. 
From the theoretical point of view, the choice of this normalisation comes from assuming 
that, in the presence of strong noise, the measure becomes smooth in the neighbourhood 
of the periodic point and therefore it is locally not different with respect to the ternary 
shift already analysed. We believe that the EI computed with respect to this normalising 
constant seems to be the relevant one for the simulations in perturbed systems since it 
does not depend on the intensity of the noise but only on parameters that are known from 
the set up of the experiment. 

The extremal index for the two cases is presented in the Figs [4] and [5] for a = 0.014 and 
a = 0.314, respectively. In both cases it is possible to see that 6 grows exponentially and 
the noise increases as e for e > 10 -6 while for very weak noise we obtain in both cases a 
result for the EI that depends on the length of the bin. This is probably due to numerical 
uncertainties intrinsic to operating in double precision. This effect can be clearly recognised 
by repeating the same experiment in single precision - results not shown here - where 9 
saturates for values of e < 10~ 4 because the effect of the single precision numerical noise 
becomes dominant. 



7. Henon attractor 

Finally, we investigate the impact of a stochastic perturbation applied onto a prototypical 
case of a map possessing a singular physical measure supported on an attractor whose 
dimension is smaller than the dimension of the phase space, d. For such dynamical systems 
in [29] we have shown numerically the convergence of maxima for the observables gi to the 
three classical EVLs whose parameters (k, <t, /i) depend on the scaling exponent of the 
measure of a ball centred around the point z from which the distance is computed. Such 
a scaling exponent turned out to be the local dimension cIl(z) of the SRB-measure. For 



this reason, the equations for the asymptotic parameters presented in Eq. 3J3 for the one 
dimensional and absolutely continuous case, should be modified as follows: 

«(0i) = O k{ 92 ) = l/{ad L {z)) K{g 3 ) = -l/{ad L {z)) *{ 9l ) = l/d L {z). (7.1) 

In [29] we have also shown that the dimension of the whole attractor can be recovered by 
averaging over different points z the local dimensions d^(z) obtained inverting the relations 



in Eq. (7.1). 
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Hereby, we choose to investigate the properties of the Henon map defined as: 

x j+ i = Vj + 1 - ax) + eCj , 12 , 
Vj+i = bxj, 

where we consider the classical parameters a = 1.4, b = 0.3 Q Recall that for the usual 
Henon attractor, points outside the basin of attraction escape to infinity. Hence, when 
the system is perturbed, it is natural to expect that, after waiting enough time, the evo- 
lution of every initial condition will escape to infinity since the noise lets the orbit explore 
the whole d- dimensional phase space. Moreover, as we have already pointed out in the 
introduction, since the physical measure is indeed smoothed by the action of the noise, so 
that the resulting attractor has the same dimension d of the phase space. Therefore, the 
EVLs' parameters depending on the local dimension d^z) of the attractor should change 
value abruptly as soon as the noise is switched on. We have selected 500 different z on 
the support of the physical measure, and a single realisation has been analysed for each of 
these z. For the deterministic dynamics, such a set up has been used to reconstruct the 
attractor's dimension by averaging the local dimension di(z) derived for each considered 
z from the EVLs [25] . Although there are no rigorous results on this system, we will such 
the numerical evidences and issues emerging from the analysis of Henon to suggest how to 
operate in a general case where theoretical results are - for some reasons - unavailable. 



In the case of purely deterministic dynamics, the EVLs parameters agree with high pre- 
cision with the theoretical predictions of [29]. The numerical experiments performed with 
varying intensity of the noise forcing (see Fig. [6J suggest that, also in this case, the sig- 
nature of the deterministic dynamics is pretty resilient. We need a rather intense noise to 
obtain a detectable smearing of the measure, such that our indicators see an absolutely 
continuous with respect to Lebesgue, invariant measure. Only when p ~ 0.1 we observe the 
escaping behaviour to infinity which causes a divergence of the EVL from the deterministic 
one with unreliable parameters, whereas in all the other cases we recover quite well the 
dependence on the dimension of the usual Henon attractor. Again, when dealing with finite 
data samples the correct reconstruction of the phase space and of the physical measure 
depends on the intensity of the noise. The case of Henon also calls the attention for the 
fact that the balls of the perturbed systems keep scaling with the local dimension when 
the system is weakly perturbed, instead of scaling with the dimension of the phase space. 
This effect is explained in Collet's paper [7] and is linked to the direction in which the 
stochastic perturbation operates: we have already seen for the shift map that the addition 
of noise in chaotic one dimensional maps does not affect much the system's behaviour at 

We remind that Benedicks and Carleson [3] proved that there exists a set of positive Lebesgue measure 
S in the parameter space such that the Henon map has a strange attractor whenever o,ieS. The value 
of b is very small and the attractor lives in a small neighbourhood of the cc-axis. For those values of a and 
b, one can prove the existence of the physical measure and of a stationary measure under additive noise, 
which is supported in the basin of attraction and that converges to the physical measure in the zero noise 
limit [5]. It is still unknown whether such results could be extended to the "historical" values that we 
consider here. 
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generic points. The effect of noise is negligible for the components parallel to the unstable 
manifold, while the effect is definitely stronger along the stable manifold. Introducing a 
noise acting as a forcing only along the stable manifold would indeed create a stronger 
smoothing effect. We will try to investigate this possibility in a future publication. The 
heuristic take-home message we learn from the case of the Henon map is that the noise 
is not an easy fix for the singularity of the invariant measure in the case of dissipative 
systems, in the sense that its effect becomes noticeable only when its intensity is relatively 
larger, unless we consider extremely long time series. 

As discussed in the introduction, this is a delicate point related to the practical applicability 
of the FDT in dissipative statistical mechanical systems. While it is clear that FDT does 
not apply for systems possessing a singular invariant measure [3U [35j [27J [30] , some authors 
[TJ |2I] have advocated the practical applicability of the FDT thanks to the smoothing effect 
of noise. It is clear that, while even an arbitrarily small noise indeed smooths the invariant 
measure, one may need to accumulate an extraordinarily long statistics record in order to 
observe it, and, therefore, to be able to apply for real the FDT. Obviously, this is a matter 
worth investigating on its own. 
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Figure 1. GEV parameters VS intensity of the noise e = 10~ p for the 
circle rotations perturbed map. Blue: n = 10 3 , m = 10 3 , Green: n = 10 4 , 
m = 10 3 . Red lines: expected values, z ~ 0.7371. From the top to the 
bottom: K,(g 2 ), K(g 3 ), a{gx). 
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Figure 2. GEV parameters VS intensity of the noise e = 10~ p for the 
ternary shift perturbed map. Blue: n = 10 3 , m = 10 3 , Green: n = 10 4 , 
m = 10 3 . Red lines: expected values, z ~ 0.7371. From the top to the 
bottom: K,(g 2 ), K,(g 3 ), a{gx). 
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Figure 3. Extremal index 9 VS intensity of the noise e = 1CU P for 
ternary shift perturbed map. Blue: n = 10 3 , m = 10 3 , Green: n = 
m = 10 3 . Red line: theoretical 6 for z = 0.5 of the unperturbed map. 
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Figure 4. Extremal index 9 (in log scale) VS intensity of the noise e = 10 _p 
for the quadratic perturbed map with a=0.314, z ~ 0.72. Blue: n = 10 3 , 
m = 10 3 , Green: n = 10 4 , m = 10 3 . 




Figure 5. Extremal index 9 (in log scale) VS intensity of the noise e = 10 _p 
for the quadratic perturbed map with a=0.014, z ~ 0.98. Blue: n = 10 3 , 
m = 10 3 , Green: n = 10 4 , m = 10 3 . 
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Figure 6. GEV parameters VS intensity of the noise e = 10~ p for the 
Henon perturbed map. Blue: n = 10 3 , m = 10 3 , Green: n = 10 4 , m = 10 3 . 
Red lines: expected values, z is different for each realisation. From the top 
to the bottom: n(g 2 ), K,{g 3 ), <r(^i). 
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